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ABSTRACT 


A literature search was conducted to examine current crack 
initiation and propagation models for notches. The Dowling crack 
initiation model and the El Haddad et al. crack propagation model 
were cho,sen for additional study. Existing data was used to make 
a preliminary evaluation of the crack propagation model. The 
results indicate that fo'^' the crack sizes in the test# the 
elastic parameter K gave good correlation for the crack growth rate 
data. Additional testing, directed specifically toward the 
problem of small cracks initiating and propagating from notches 
is necessary to make a full evaluation of these initiation and 
propagation models. 
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FOREWORD 


This First Interim Technical Report is submitted in 
accordance with the NASA Grant MAG ’3-246. It summarizes the 
work accomplished during the period 1 January 1982 through 
31 August 1982. 

The work is being performed for the National Aeronautics 
and Space Administration. The work is being performed by the 
Aerospace Mechanics Department of the University of Dayton 
Research Institute, with Mr. George Hartman and Dr. J. P. Gallagher 
acting as co-Principal Investigator. Technical support is being 
provided by Dr. A. M. Rajendran and Mr. D. S. Dawicke. 
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SECTION 1 
INTRODUCTION 


The growth of cracks can be divided into the two categories 
of initiation and propagation. In a notched structure, the 
initiation and subsequently the propagation of small cracks are 
governed by the stress field created by the notch. After the 
crack grows beyond the influence of the notcn, crack propagation 
behavior is governed by the nominal stress in the structure. 

A literature search was performed to examine the state-of- 
the-art crack propagation and initiation models. Thirteen recently 
published studies [1-13] provide examples of crack initiation and 
propagation models. The most promising of these models are: 
the El Haddad et al. crack propagation model [1], the elastic- 
plastic J-integral based crack propagation model described by 
Ratwani et al. [2], and the Dowling crack initiation cut-off 
model based on linear elastic fracture mechanics [3,4,5] . The 
work to date has included a close examination of the promising 
crack propagation models, specifically by testing their ability 
to describe existing data. 


SECTION 2 
LITERATURE REVIEW 


A literature review was conducted to determine the 
existing, approaches taken to solve the problem of cracks initi- 
ating and propagating from notches. This re^'iew was presented 
at NASA-Lewis on 16 July 1982, The remainder of this section 
summariaes the important /Aspects of the crack initiation an-'l 
propagation models that have been chosen for additional study. 

2.1 STRAIN LIFE INITIATION MODEL 

A crack initiation model developed by Dowling [3,4,5} 
predicts life to initiation for notch geometries using smooth 
specimen strain life data. The severity of ^the notch is taken 
into account by elevating the applied stress by the elastic 
stress concentration factor (k^) . The equivalent strain ampli- 
tude for notched geometry is obtained from the cyclic stress- 
strain curve of the material. 

Dowling found this method to be effective in predicting 
crack initiation life if the initiation crack length remained 
within the influence of the notch. Dowling defined an upper 
bound on the initiation crack length which ensured the crack was 
within the influence of the notch. This upper bound was defined 
to be the intersection of what he termed the short crack stress- 
intensity factor (Kg) and the long crack stress-intensity factor 
(K^) solutions described as a function of crack length. The 
short crack stress-intensity factor was defined as a function of 
the stress concentration factor (k^) of the notch and the crack 
length measured from the edge of the notch (il) , 

K = 1.12S k, /il (1) 

s t 

where 1.12 is the free surface correction factor for an edge 
crack. The long crack stress-intensity factor assumed the notch 




has no effect on the crack behavior, 

w S /ia (2) 

where a is the sum of the crack length' (Jt) and radius (r) of the 
notch; s in Equations 1 and 2 is the applied stress. 

long crack and short crack stress- intensity factor 
solutions for a diametrically cracked 1.0 inch diameter hole in a 
wide plate are shovm in Figure 1. Analytically, the crack length 
associated with the intersection of the long crack and short 
crack stress-intensity factors can be expressed in terms of notch 
radius (r) and stress concentration factor of the notch. 


(1.12 k^)^-l 

The Dowling crack initiation cut-off parameter given by 

Equation 3 is shown in Figure 2 as a function of notch severity. 

As the notch severity increases the initiation crack length 

# 

decreases, i.e., the notch behaves like a crack. It can also be 
noted that the initiation crack length is on the order of 10 to 
20% of the notch radius. 

Dowling also suggested that linear elastic fracture mecha- 
nics parameters could be utilized beyond the initiation cut-off. 

2.2 CRACK PROPAGATION MODELS 

2.2.1 Elastic J-Integral (K) 

The linear elastic fracture mechanics (LEFM) 
parameter K, the stress-intensity factor, has been found to be 
an effective basis for life analysis procedures when extensive 
plasticity does not occur. Localized cyclic plasticity around 
the root of a notch has been observed to accelerate the growth 
rate of cracks propagating from notches. El Haddad et al. [1], 
Leis et al. [11,12], and Newman [13] have shown that the LEFM 
parameter K, correlates the crack growth rate data in a manner 
which would lead to unconservative life predictions for small 
cracks propagating from plastically deforming notches under fully 
reversed loading. 
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If __ 



i (in.) 

Crack length measured from edge of notch. 


Figure 1. The Long and Short Crack Stress-Intensity Factor 
Solutions for a 1 inch Diametrically Cracked Hole. 






Dowling’s crack initiation cutoff parameter » 


Figure 2. Effect of Notch Severity on Dowling’s Crack 
Initiation Cutoff Parameter. 



The advantage of using the LEFM parameter K, is 
that it is a material independent property of a specific geometry. 
Handbook solutions exist for many simple geometries and a series 
of relatively inexpensive linear finite element runs will generate 
the solution for more complex geometries. One of the more direct 
methods ior calculating the stress-intensity factor via finite 
element analysis involves using the path independent line integral 
derived by Rice [14]^ i.e., the J-integral, which is given by 

J - / (Wdy - T || dB) (4) 

In the elastic case, the stress-intensity factor is proportional 
to the square root of the J-integral according to the formula, 

K « /JK (5) 

where E is the elastic modulus. The LEFM parameter K provides a 
starting point in the life analysis of cracks propagating from 
notches due to the availability of the solution and its acceptance 
-In making linear elastic life predictions. 

2.2.2 Elastic-Plastic J-Integral 

Ratwani et al. [2] suggested that the behavior of 
cracks propagating from plastically deforming notches could be 
described using the j-integral calculated to account for the 
plastic behavior. The J-integral proposed for this correlation 
is calculated with Equation 4 based on finite element results for 
elastic-plastic material properties. Actually, once Ratwani et 
al. calculated the value of the J-integral under elastic-plastic 
conditions for a given loading and crack length, they converted 

J to a quasi-stress intensity factor (K ) using Equation 5, i.e., 

It 

K = ✓JB (6)_ 

which was then used for calculating crack growth behavior in the 
notch region. 
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2.2.3 


Strain Based Stress-.Tntensity Factor 


El Haddad et al. [1] proposed a different type of 

a ' 

plasticity enhanced stress-intensity factor to account for the 
behavior of cracks propagating from notches. This approach is 
based upon a strain concentration factor (k' ) , which is the ratio 
of loca3^- strain (e) as a function of distance (x) from the notch 
to the nominal gross strain (e) . 

kj. = (X) = (7) 

The stress-intensity factor (K^j,) proposed by El Haddad et al. is 
given by 

Kgj, = s (8) 

where £ is the crack length measured from the notch root and k^ 
is evaluated for x«A (See Equation 7) . 

The strain concentration factor is evaluated by 
a nonlinear finite element solution of the uncracked* geometry . 

The El Haddad et al. stress-intensity factor is equivalent to the 
LEFM K for different levels of load when the material at the notch 
root does not experience plastic deformation. .For stress levels 
which result in plastic deformation at the notch root, the strain 
concentration factor will elevate the stress-intensity factor to 
account for the accelerated growth rate of the crack. For crack 
lengths outside the plastic zone of the notch, the stress-intensity 
factor approaches that of the LEFM long crack stress-intensity 
factor. El Haddad et al. presented results which indicate this 
method correlates the crack growth rate data for cracks propa- 
gating from several types of notches under several loading 
conditions , 
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SECTION 3 

APPLICATION OF EXISTING DATA TO 
CRACK PROPAGATION MODELS 


A small data base of crack growth rate data [15] and a 

series of finite element results [16] for industrial grade buss 

* 

bar copper were available for use in analyzing the crack propa- 
gation models. Fatigue crack growth rate data had been collected 
using 4.0 inch wide plates with center cracks (CCP) and with 
radially cracked holes (RCH) . Specimen types and test stress 
levels are given in Table 1. Tensile tests and incremental step 
tests were conducted to obtain the monotonic and cyclic stress- 
strain curves for material modeling. The finite element based 
study was conducted using the MAGNA code and the material model was 
obtained from the tensile test results. 


3.1 ELASTIC J-INTEGRAL (K) 

A series of elastic finite element solutions wfere obtained 
for a crack propagating from the RCH specimen. The LEFM para- 
meter K, stress-intensity factor, for this geometry was calculated 
from the J-integral results of the elastic finite element solu- 
tion. Figure 3 describes the finite element grids, crack tip 
elements, and a J-integral path for a 1 inch crack as measured 
from the edge of the 1 inch diameter, hole in the 4 „inch wide plate. 
The finite element solution for the stress-intensity factor is 
shown in Figure 4. The Bowie [17] solution for the stress- 
intensity factor of a crack propagating from a 1 inch diameter 
hole in an infinite plate has been added for comparison. 

The difference between the Bowie solution and the finite 
element solution for K is due to the finite width of the plate. 

The RCH crack growth rate data correlated with the finite 
element determined K are shown in Figure 5, along with the crack , 
growth rate data for the CCP specimens. The fatigue grovjth 

rate obtained from these specimens were correlated with the 
analytical elastic K solution with a width correction factor used 
in A^TM standard E647 [18] . 
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CRACK GROWTH RATE TESTS FOR INDUSTRIAL 
GR7UDE BAR COPPER TESTED AT A STRESS 
RATIO (R) OF 0.1 
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Coarse Finite Element Mesh 



Fine Finite Element Mesh for Crack Tip 
Region 


Figure 3. Finite Element Mesh Representation of a 4 inch 
Wide Plate with a Crack Propagating from One 
Side of a 1 inch Diameter Hole. 
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Figure 4 . The Elastic Finite Element Solution and the 
Bowie Approximation. 
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Figure 5. Crack Growth Rate Data Correlated w 
Elastic Stress-Intensity Factor for 
Industrial Grade Buss Bar Copper Su 
to a Fatigue Loading with Stress Ra 
of 0.1. 
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(9) 


K ® S ^Tia Sec 

where a is the. half crack length and W is the plate width. 

For this material, the correlation with the finite element 
elastic K is very good. The accelerated crack growth rate for 
small cracks near the notch reported by others was not observed. 
The smallest crack len<'^ths shown in Figure 5 are on the order of 
0.1 inch which is small compared to the calculated size of the 
plastic zone shown in Figure 6. 

A study by Newirtan [131 suggests that crack closure is the 
cause for the accelerated crack growth rate experienced by small 
cracks propagating from notches. The current data base seems to 
agree with this conclusion. The previously reported accelerated 
crack growth rate was seen in tests run with a stress ratio of 
-1 (fully reversed loading) . The existing test results were for 
tests at a stress ratio of 0,1, thus crack closure may not be as 
important in the cracking process. We are planning further 
investigation in this area. 

Newman reasoned that once closure has occurred, the load 
needed to open a short crack is less than that needed to open a 
long crack. This would result in a greater stress range in which 
damage is being done for the short crack, thus an increased crack 
growth rate. 

3.2 ELASTIC-PLASTIC J-INTEGRPi 

Nonlinear finite element based studies were also conducted 
for cracks propagating in the RCH specimens. Figure 7 summarizes 
the results of the J-integral calculations in terms of the ratio 
of "plastically" corrected stress-intensity factor K (given by 

IT 

Equation 6) to the elastic stress-intensity factor, for the two “ 
applied stress levels utilized during the test. 
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Plastic Zone Size As Hsasured From Notch (In.) 


Figure 6. Extent of Plastic Zone Along Expected Path of 
Crack for a 1 inch L/iameter Hole in a 4 inch 
Wide Plate. 
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Figure 7. Extent of Plastic Enhancement of the Stress- 
Intensity Factor, as Calculated by the Finite 
Element Solution. 



The '’plastically” corrected stress-intensity factors based 
on Equation 6 were used to correlate the crack growth rate data 
from the RCH specimens as shown in Figure 8, Finite element 
solutions for the CCP specimens are not complete; thus, the crack 
growth rate data for the CCP specimens were not included in the 
figure. { 

$ 

3.3 STRAIN BASED STRESS- INTENSITY FACTOR 

A nonlinear finite element solution for an uncracked RCH 
specimen was used as an initial estimate of the plastic strain 
near the notch. The results of this program were used to calcu- 
late the strain concentration factor (ki) . This factor was used 

e 

in the plastically enhanced stress-intensity parameter . 

KgH then used to correlate the crack growth rate data for the 
four RCH tests. The results of this correlation, along with the 
elastic long crack (CCP) data are shown in Figure 9. 

i 

The plastically enhanced stress-intensity fact'or, as 
described by El Haddad et al., failed to correlate the crack 
growth rate data from the RCH specimens with the long crack data. 
An additional test was conducted to determine the validity of 
the strains calculated by the finite element analysis. The test 
specimen consisted of an uncracked RCH specimen with strain gages 
attached to the inside of the notch. The specimen was then cycled 
at a stress level similar to that used in the fatigue tests. The 
experimental results indicate that the material at the notch root 
experienced a large decrease in plastic strain after the first 
cycle due to redistribution of stresses. The plastic strain 
amplitude stabilized after approximately the tenth cycle. The 
finite element solution agrees with the first cycle but differs 
substantially with the tenth and succeeding cycles, as shown in - 
Figure 10. The strain concentration factor for the uncracked 
RCH specimen, as calculated by the finite element solution, does 
not give an accurate description of the strain along the expected 
crack path. A cyclic nonlinear finite element analysis is 
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Figure 8. Crack Growth Rate Data Cor: 
Enhanced Stress-Intensity ] 
r Solution. 
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Figure 10. Fxtent of Redistribution of Stresses Experienced 
by the Copper and Finite Bl^jjnent Predictions of 
Strain. 





planned and when the data is available it will be u«ed to 
obtain a better estimate o£ the plastic strain for the El Haddard 

a 

model * 

3.4 SUMMARY OF EXPERIMENTAL RESULTS 

The experimental results used thus far to evaluate the 
crack propagation models show that the elastic stress^^intensity 
factor (K) correlates the crack growth rate data for the condi- 
tions tested. The elastic-plastic J-integral based (Kp) and the El 
Haddad et al, plastically enhanced stress-intensity factor 
are strongly dependent upon the material model used in the finite 
element analysis. 

Additional testing is 4.equired to determine the conditions 
under which the elastic K fails to correlate the data. Tests are 
planned which will specifically address the problem of cracks 
propagating from notch, es. These tests will provide crack growth 
rate data for small cracks (less than 0.05 inch) near notches. 

Tests will be run under fully reversed loading and under a stress 
ratio of 0.1, to determine the effect of crack closure on the 
small crack growth rate. 


SECTION 4 


SUMMARY OP WORK ACCOMPLISHED AND WORK 
PLANNED FOR THE NEXT PERIOD 

During the past six months, the work accomplished can be 
summarized as follows: 

A. ' A literature review of current crack initiation and 
propagation model for notched structures was conducted. 

B. A crack initiation model and several crack propagation 
models were chosen for further study. 

C. Analysis, using existing crack growth rate data and 
finite element solutions, was performed to evaluate the crack 
propagation models. 

D. The materials for additional testing has been acquired 
and the machining of the specimens has begun. 

Between September 1 and 31 December 1982, the following 
work will be accomplished: 

A. Baseline testing of the aluminum will be conducted. 

B. Life predictions using the crack initiation and ctack 
propagation models will be made. 

C. Fatigue tests will be conducted on the aluminum. 

D. The crack initiation and crack propagation models will 
be evaluated for the aluminum. 

E. A plate of INCONEL 718 will be procured. 
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